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Thermoelectric  systems  (TE)  can  directly  convert  heat  to  electricity  and  vice-versa  by  using  semiconductor  ma¬ 
terials.  Therefore,  coupling  between  heat  transfer  and  electric  field  potential  is  important  to  predict  the  perfor¬ 
mance  of  thermoelectric  generator  (TEG)  systems.  This  paper  develops  a  general  two-dimensional  numerical 
model  of  a  TEG  system  using  nanostructured  thermoelectric  semiconductor  materials.  A  TEG  with  p-type  nano¬ 
structured  material  of  Bismuth  Antimony  Telluride  ( BiSbTe )  and  n-type  Bismuth  Telluride  ( Bi2Te3 )  with  0.1  vol.% 
Silicon  Carbide  (SiC)  nanoparticles  is  considered  for  performance  evaluations.  Coupled  TE  equations  with  tem¬ 
perature  dependant  transport  properties  are  used  after  incorporating  Fourier  heat  conduction,  Joule  heating, 
Seebeck  effect,  Peltier  effect,  and  Thomson  effect.  The  effects  of  temperature  difference  between  the  hot  and 
cold  junctions  and  surface  to  surrounding  convective  on  different  output  parameters  (e.g„  thermal  and  electric 
fields,  power  generation,  thermal  efficiency,  and  current)  are  studied.  Selected  results  obtained  from  current  nu¬ 
merical  analysis  are  compared  with  the  results  obtained  from  analytical  model  available  in  the  literature.  There  is 
a  good  agreement  between  the  numerical  and  analytical  results.  The  numerical  results  show  that  as  temperature 
difference  increases  output  power  and  amount  of  current  generated  increase.  Moreover,  it  is  quite  apparent  that 
convective  boundary  condition  deteriorates  the  performance  of  TEG. 
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1.  Introduction 

The  research  and  development  in  nanostructured  TE  systems  have 
gathered  considerable  attention  due  to  their  potential  applications  in  di¬ 
rect  electricity  generation,  refrigeration,  and  air-conditioning.  TE  systems 
can  largely  be  classified  as  thermoelectric  generator  (TEG)  and  thermo¬ 
electric  cooler  (TEC).  The  TEG  converts  heat  into  electricity  and  TEC  con¬ 
verts  electricity  into  heating/cooling  based  on  Seebeck  and  Peltier  effects, 
respectively.  TE  systems  are  solid  state  heat  engines/refrigerators  which 
are  robust,  silent,  compact,  and  environment  friendly.  TE  systems  are 
made  up  of  numbers  of  p-type  and  n-type  semiconductor  elements  con¬ 
nected  electrically  in  series  and  thermally  in  parallel.  The  TEC  has  a 
wide  range  of  applications;  for  example,  electronic  cooling,  laser  diode 
cooling,  military  garment,  laboratory  cold  plates,  and  automobile  seat 
cooler.  In  a  similar  manner,  TEG  has  various  applications  in  military, 
deep  space  vehicles,  remote  power  sources  for  inhabitable  places,  solar 
and  waste  heat  power  generator.  Liquid  cooling  of  CPU  using  TE  was  pro¬ 
posed  and  experimentally  investigated  considering  different  material  and 
size  of  the  heat  sinks  [1  ].  Huang  et  al.  [2]  discussed  TE  cooling  of  electronic 
equipment  experimentally  and  analytically.  The  results  determined  that 
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the  integration  of  water-cooling  with  TE  is  helpful  to  increase  the  perfor¬ 
mance  of  electronic  equipment.  In  order  to  address  the  site  specific  on- 
demand  cooling  of  hot  spot  in  microprocessor  [3],  a  numerical  simulation 
that  includes  heat  spreader,  thermal  interface  material,  chip,  and  nine 
TECs  was  carried  out  considering  steady  state  and  transient  analysis.  Sul¬ 
livan  et  al.  [3]  concluded  that  transient  cooling  with  square  root  current 
pulse  is  most  effective  with  10  °C  cooling.  Wang  [4]  proposed  and  inves¬ 
tigated  experimentally  the  TEG  using  waste  heat  of  the  Light-Emitting  Di¬ 
odes.  Results  reported  by  Wang  [4]  investigated  power  output  of  TEG 
using  waste  heat  of  the  Light-Emitting  Diodes  (LED).  Experiments  of 
Wang  [4]  showed  160  mW  of  power  output  from  TEG  with  6  W  of 
input  power  to  LED.  Recently,  Hsiao  et  al.  [5]  studied  the  performance  of 
TE  modules  as  a  waste  heat  recovery  tool  from  an  automobile  engine 
using  a  one-dimensional  thermal  resistance  model  and  compared  their 
model  with  experimental  data.  Results  [5]  showed  that  the  performance 
of  a  TE  module  on  the  exhaust  pipe  performs  better  compared  to  a  TE 
module  on  the  radiator  system.  Rezania  et  al.  [6]  studied  the  effect  of 
cooling  power  on  the  performance  of  a  TEG.  Rezania  et  al.  [6]  determined 
the  optimum  coolant  flow  rate  for  maximum  power  output  for  the  TEG. 
For  example,  temperature  difference  of  10  K  gives  maximum  power  out¬ 
put  of  0.035  W  with  coolant  flow  rate  of  0.07  1/min  [6],  One-dimensional 
analytical  solutions  of  conventional,  composite,  and  integrated  TEG  were 
carried  out  by  Reddy  et  al.  [7]  considering  adiabatic  and  convective  side 
wall  conditions.  Regardless  of  TEG  design,  an  increment  in  the  hot  side 
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Nomenclature 

C  specific  heat  capacity  (kJ/kgK) 

D  electric  flux  density  vector  (C/m2) 

E  electric  field  intensity  vector  (N/C) 

e  electron  charge  (1.60  x  10_19C) 

h  convection  heat  transfer  coefficient  (W/m2I<) 

H  height  of  TE  leg  (m) 

I  current  (A) 

J  electric  current  density  (A/m2) 

k  Boltzmann  constant  (JK-1 ) 

L  length  (m) 

L„  Lorentz  number  (V2K-2) 

R  resistance  (fi) 

q  heat  (W) 

q  heat  generation  rate  per  unit  volume  (W/m3) 

Q  heat  (W) 

t  time  (S) 

T  temperature  (K) 

W  width  of  the  TE  leg  (m) 

x  coordinate  (m) 

y  coordinate  (m) 

ZT  dimensionless  figure  of  merit 

p  p-type  semiconductor  material 

n  n-type  semiconductor  material 


Greek  symbols 

a  Seebeck  coefficient  (V/K) 

p  density  (kg/m3) 

A  thermal  conductivity  (W/mK) 

a  electrical  conductivity  (S/m) 

4>  electric  scalar  potential  (V) 

n  thermal  efficiency 

e  dielectric  permittivity  matrix  (F/m) 


Subscripts 

C  cold  surface 

corn  convection 

E  electronic 

g  gap 

H  hot  surface 

L  external  load 

L  lattice 

p  p-type  material 

n  n-type  material 


number  and  fluid  temperature  on  performance  of  an  integrated  TEG  sys¬ 
tem.  Reddy  et  al.  [10]  found  that  higher  Reynolds  number  enhances  heat 
transfer  and  thus  leads  to  higher  power  output  of  TEG.  Zhou  et  al.  [1 1  ]  de¬ 
veloped  simple  and  coupled  field  model  with  the  former  considering 
Navier-Stokes  and  energy  equations  with  continuity  equation  and  latter 
with  different  TE  effects  such  as  Seebeck,  Peltier,  and  Thomson  effect. 
Zhou  et  al.  [1 1  ]  reported  overall  TEG  efficiency  of  3.5%  with  temperature 
difference  of  80  K.  Baranowski  et  al.  [12]  developed  mathematical 
model  for  solar  TEG  which  can  provide  analytical  solutions  of  device  effi¬ 
ciency  with  temperature  dependant  properties.  They  [12]  also  showed 
that  considering  currently  available  materials,  total  efficiency  of  14.1%  is 
possible  with  cold  and  hot  side  temperatures  of  100  °C  and  1000  °C,  re¬ 
spectively.  Baranowski  et  al.  [12]  also  concluded  that  if  figure  of  merit 
(ZT)  reaches  2  then  TEG  efficiency  of  25%  can  be  obtained,  for  the  cold 
and  hot  side  temperatures  of  100  °C  and  1000  °C,  respectively. 

Altenkirch  [13]  introduced  figure  of  merit,  ZT  as  a  parameter  to  clas¬ 
sify  different  materials.  The  performance  of  TE  materials  are  character¬ 
ized  by  a  dimensionless  parameter  ‘figure  of  merit',  ZT  =  ( o?a/K)T , 
where  a  is  the  Seebeck  coefficient,  cris  the  electrical  conductivity,  and 
\  is  the  thermal  conductivity.  The  current  ZT  value  of  the  best  available 
TE  materials  is  1  at  room  temperature  [14],  Slack  [15]  described  that 
good  TE  materials  need  to  have  low  thermal  conductivity.  In  addition 
to  this,  Slack  [15]  described  that  the  best  TE  material  would  behave  as 
“Phonon  Glass  Electron  Crystal  (PGEC)";  that  is,  it  would  have  thermal 
properties  of  glass  like  material  and  electrical  properties  of  crystalline 
material.  Recent  advancements  in  the  field  of  nanotechnology  [14] 
have  opened  the  door  for  further  improvements  of  ZT  for  the  TE  mate¬ 
rials.  The  expression  of  the  figure  of  merit,  ZT  =  (a^cr/X)  T,  evidently  in¬ 
dicates  that  one  of  the  methods  to  increase  figure  of  merit  is  to  reduce 
the  thermal  conductivity  of  the  TE  material.  Thermal  conductivity  is 
the  sum  of  two  contributions:  electrons  and  hole  transporting  the  heat 
(X£)  and  phonons  traveling  through  the  lattice  (X£)  [16],  The  electronic 
part  of  thermal  conductivity  (X£)  is  related  to  the  electrical  conductivity 
as  per  the  Wiedemann-Franz  law  [16]  as  shown  in  Eq.  (1), 


where,  La  is  the  Lorentz  number  and  for  metals  it  is  equal  to  [16], 

n2k2  s  ,  •, 

T0  =  “^5-  =  2.45  x  10-8V2K-2  (2) 

where,  k  is  the  Boltzmann  constant  ( 1 . 3  8  x  1 0  “  23  JK_  1 )  [  1 6  ]  and  e  is  the 
electron  charge  (1.60  x  10_19C)  [17],  The  expression  of  the  figure  of 
merit  can  be  written  in  terms  of  Lattice  conductivity  (\L)  and  electronic 
conductivity  (X£)  as  shown  below  [16], 


ZT- 


a2  X£ 

Loh+K 


(3) 


temperature  enhances  the  performance  of  TEG  systems  [7],  Reddy  et  al. 
[7]  concluded  that  composites  and  integrated  TEG  extracts  more  heat 
compared  to  conventional  TEG  and  reduces  rare-element  material 
usage.  Gou  et  al.  [8]  investigated  the  performance  of  low  temperature 
waste  heat  TEG  using  one-dimensional  analytical  simulations  and  exper¬ 
iments.  They  [8]  concluded  that  in  addition  to  increasing  the  waste  heat 
temperature  and  number  of  modules  in  series,  expanding  the  heat  sink 
surface  area  and  enhancing  the  cold  side  heat  transfer  in  proper  ranges 
can  have  dramatic  effects  on  TEG’s  performance.  A  three  dimensional  nu¬ 
merical  model  of  TEG  applied  to  fluid  power  systems  was  developed  by 
Chen  et  al.  [9]  and  their  numerical  simulation  performed  with  ANSYS 
showed  fairly  good  match  with  experimental  results.  In  addition  to  this, 
Chen  et  al.  [9]  concluded  that  convection  heat  transfer  losses  increase 
the  heat  input  to  TEG  thus  reducing  the  thermal  efficiency  of  TEG.  A 
three-dimensional  coupled  numerical  simulation  of  integrated  TE  device 
was  carried  out  by  Reddy  et  al.  [10]  to  check  the  effects  of  Reynolds 


One  method  to  enhance  ZT  is  the  inclusion  of  nanoparticles  into  the 
bulk  TE  materials  which  can  lead  to  low  lattice  thermal  conductivity 
[18],  Poudel  et  al.  [19]  achieved  ZT  value  of  1.4  at  373 K  by  hot-pressing 
of  nanopowders  of  Bi2Te3  and  Sb2Te3  under  Argon  (Ar)  atmosphere.  The 
enhanced  ZT  was  due  to  the  significant  decrease  in  the  lattice  thermal 
conductivity  of  material.  Zhao  et  al.  [20]  fabricated  Bi2Te3  with  various 
amounts  of  nano  SiC  particles  using  mechanical  alloying  and  spark  plasma 
sintering  and  tested  the  TE  and  mechanical  properties.  Their  results 
showed  an  improvement  in  ZT  from  0.99  to  1.04  with  an  inclusion  of 
0.1%  vol.%  SiC  particles.  Li  et  al.  [21]  obtained  ZT  of  1.43  for  double- 
doped  C04Sbi2  skutturedites  using  Indium  (In)  and  Cerium  (Ce)  doping. 
The  attractive  results  achieved  using  nanotechnology  has  encouraged  re¬ 
searchers  to  include  nanopartide-doped  TE  materials  for  various  low  po¬ 
tential  heat  recovery  applications;  for  example,  solar  TEG  and  automobile 
exhaust  heat  recovery.  Kraemer  et  al.  [22]  have  proposed  novel  solar  TEG 
with  glass  vacuum  enclosure  considering  nanostructured  TE  materials. 
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The  developed  solar  TEG  achieved  maximum  efficiency  of  4.6%  with  solar 
flux  of  1000  W/m2  condition.  McEnaney  et  al.  [23]  developed  a  novel  TEG. 
They  [23]  placed  high  performance  nanostructured  material  in  evacuated 
tube  with  selective  absorber  and  achieved  an  efficiency  of  5.2%.  It  is  quite 
evident  that  nanostructured  TE  materials  can  increase  the  performance  of 
TE  systems. 

It  is  quite  clear  from  the  above  discussion  that  in  future  nanostruc¬ 
tured  TE  materials  will  play  a  significant  role  as  a  direct  energy  conversion 
tool  from  low  potential  sources.  It  is  necessary  to  understand  the  wide 
range  of  characteristic  features  of  newly  developed  TE  materials.  Such 
characteristic  features  include  heat  transport  and  electric  potential. 
The  existing  literature  on  the  numerical  simulation  of  nanostructured 
TEG  is  very  limited.  In  this  paper,  performance  of  nanostructured  TEG  is 
evaluated  using  2D  numerical  simulation.  The  temperature  dependent 
thermophysical  and  electrical  properties  of  nanostructured  TE  material, 
surface  to  surrounding  convection  heat  transfer  losses,  and  Thomson  ef¬ 
fect  are  included  in  the  current  model.  The  field  plots  of  heat  and  current 
are  presented  with  different  convection  heat  transfer  coefficients.  Numer¬ 
ical  results  are  compared  with  that  of  one  dimensional  analytical  results 
in  terms  of  current  produced. 


heat,  temperature,  heat  generation  rate  per  unit  volume,  and  heat  flux 
vector,  respectively.  The  continuity  of  electric  charge  through  the  sys¬ 
tem  must  satisfy, 


(2) 


where  J  is  the  electric  current  density  vector  and  D  is  the  electric  flux 
density  vector,  respectively.  Eqs.  (1)  and  (2)  are  coupled  by  the  set  of 
TE  constitutive  equations  [24]  as  shown  in  Eqs.  (3)  and  (4), 

q  =  T[a].J-[\].VT  (3) 


J  =  [O]  -  (E— [aj  -  V7-)  (4) 

where  [a]  is  the  Seebeck  coefficient  matrix,  [X]  is  the  thermal  conductiv¬ 
ity  matrix,  [a]  is  the  electrical  conductivity  matrix,  and  E  is  the  electric 
field  intensity  vector,  respectively.  E  can  be  expressed  as  —  V</>,  where 
( i >  is  the  electric  scalar  potential  [25],  Rearrangement  of  Eqs.  (1)  to  (4) 
gives  the  coupled  TE  equations  for  heat  transfer  and  electric  potential  as, 


2.  Mathematical  model  and  boundary  conditions 

The  two-dimensional  schematic  diagram  of  the  TEG  being  investi¬ 
gated  is  shown  in  Fig.  1.  The  TEG  is  mainly  comprised  of  two  vertical 
p-type  and  n-type  semiconductor  legs  connected  electrically  in  series 
and  thermally  in  parallel.  Both  legs  are  connected  through  an  electrical¬ 
ly  conductive  copper  strip.  Each  leg  has  a  cross-sectional  area  of  LxW, 
height  of  H,  and  separated  by  distance  Lg  as  shown  in  Fig.  1.  The  Qh  and 
Qc  are  amount  of  heat  available  at  heat  source  and  heat  sink  respective¬ 
ly.  The  Qjconv  is  convection  heat  loss  through  the  side  walls  of  the  TEG. 
During  the  analysis  following  assumptions  are  made: 

•  TE  materials  are  homogeneous  and  isotropic  with  temperature  de¬ 
pendent  properties. 

•  Contact  resistances  at  interface  of  copper  and  TE  materials  are 
neglected. 

In  the  TEG,  the  energy  transport  and  current  flow  are  governed  by 
energy  equation  and  continuity  of  current  density  as  per  below, 

pc|  +  V.q=%m  (1) 


pC"i  +  V  •  (T[a]  •  J)  -  V  -  ([\]  •  VT)  =  E  ■  j  (5) 

and 

V '  ([s] '  w) +  v ' ( [<7] ' W  ■ ' vr)  +  V ' ( [o] ' v<w  =  0  (6) 

where  [e]  is  the  dielectric  permittivity  matrix  [25]  and  E  J  represents 
Joule  heat  [24], 

In  Eq.  (5),  the  second  term  represents  the  Thomson  heat  and  the 
third  term  represents  heat  transfer  due  to  conduction.  In  Eq.  (6)  the 
first  term  represents  electric  current  density  due  to  the  Seebeck  effect 
and  standard  voltage  driven  electric  current. 

Thermal  boundary  conditions  for  (Fig.  1)  are  as  follows; 

•  The  top  surface  of  the  TEG  experiences  constant  hot  temperature  (TH). 

•  The  bottom  surface  experiences  constant  cold  temperature  (Tc). 

•  The  vertical  surfaces  of  p  —  type  and  n  —  type  are  considered  with  two 
different  conditions:  convective  heat  transfer  condition  and  adiabatic 
condition  (special  case). 


where  symbols  p,  Cp,  T,  q,  and(£en  represent  material  density,  specific 


Fig. 


1.  Schematic  of  unit  cell  of  TEG. 


3.  Results  and  discussion 

In  this  section,  the  performance  of  the  unit  cell  of  TEG  is  investigated 
based  on  the  results  obtained  by  numerically  solving  governing  equa¬ 
tions  presented  in  the  previous  section.  For  p-type  material,  the  nano¬ 
structured  semiconductor  Bismuth  Antimony  Telluride  ( BiSbTe )  is 
considered  to  analyze  the  performance  [19],  While  for  n-type  material, 

Table  1 

Polynomial  functions  of  Seebeck  coefficient,  electrical  conductivity,  and  thermal 
conductivity  as  a  function  of  temperature  for  BiSbTe  nanostructured  bulk  alloys  and 
Bi2Te3  with  SiC  nanoparticles  [19,20], 


Properties  Polynomial  expressions 

-139.41  +  1.43T  -  5.71  x  10^3T2  +  10.71  x  10_ST3  -  6.80  x  10 
-*!* 

-338.53  +  7.13T  -  2.73  x  10  2T2  +  4.15  x  10  T5  -  2.23  x  10 

-T* 

5.82  -  0.0323T  +  1.035  x  IQ-4!2  -  1.61  x  10  7P  +  9.64  x  10 
“I'T4 

-4.98  x  103  +  47.8T  -  1.73  x  10  *1*  +  2.73  x  10  4T!  -  1.59  x  10 
7 14 

9.79  -  8.20  x  10  2T  +  3.03  x  10  4T2  -  5.05  x  10  7T3  +  3.11  x  10 
~10T4 

8.25  -  653  x  10  2T  +  2.16  x  10  4T2  -  3.22  x  107T3  +  1.85  x  10 
-w14 
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Fig.  2.  (a)  Contours  of  heat  transfer  and  (b)  contours  of  electric  potential  with  adiabatic  heat  transfer  condition  (h  ss  0  W/m2K). 


Bismuth  Telluride  (Bi2Te3)  with  nano-particles  of  silicon  carbide  (SiC)  is 
considered  to  analyze  the  performance  [20],  The  temperature  depen¬ 
dent  transport  parameters  Seebeck  coefficient  (a),  electrical  conductiv¬ 
ity  (cr),  and  thermal  conductivity  (A.)  are  specified  as  polynomial 
functions  of  temperature  as  shown  in  Table  1 .  These  properties  are  eval¬ 
uated  at  an  average  temperature  of  working  temperature  range.  The 
semiconductor  leg  of  unit  cell  of  TEG  has  dimensions  of  1.5  mm 
x  1.5  mm  and  height  of  5  mm.  The  gap  between  two  consecutive 
legs  is  0.3  mm.  The  terminals  of  bothp-type  and  n-type  semiconductor 
legs  are  connected  with  external  load  (Rl)  which  are  matched  with  total 
internal  resistance  of  the  TEG.  Figs.  2(a),  3(a),  4(a),  and  5(a)  show  field 
plots  of  temperature  and  heat  flow.  Temperature  is  presented  by 
marked  isothermal  lines  with  multi-colored  background,  while  the 
heat  flow  is  presented  by  the  vertical  lines  with  arrows.  Similarly, 
Figs.  2(b),  3(b),  4(b),  and  5(b)  present  the  field  results  of  electric  scalar 
potential  and  current  flow.  The  electric  potential  is  indicated  by  the 
marked  iso-potential  lines  with  multi-colored  background,  while  the 


current  flow  is  indicated  by  the  lines  with  arrows.  Surface  to  surround¬ 
ing  convection  heat  losses  are  also  considered  from  vertical  walls  of 
both  semiconductor  legs  ( h  =  15-50  W/m2K). 

The  temperature  remains  nearly  constant  at  a  given  location  of  dis¬ 
tance  x  when  convection  is  absent  (i.e.,  h  =  0),  irrespective  of  the  p- 
type  or  n-type  semiconductor  leg.  However,  the  location  of  the  same 
isothermal  line  is  changed  in  p-type  and  n-type  legs  due  to  the  different 
transport  properties;  more  specifically,  the  thermal  conductivities.  Due 
to  the  absence  of  the  surface  to  surrounding  convection,  the  heat  flux 
lines  are  parallel  in  bothp-type  and  n-type  legs  and  uniformly  distribut¬ 
ed  over  the  cross-section  of  the  legs.  An  introduction  of  the  surface  to 
the  surrounding  convection  makes  the  temperature  distribution  non¬ 
linear  as  observed  in  Figs.  3(a),  4(a),  and  5(a).  An  increment  in  convec¬ 
tion  heat  transfer  coefficient  increases  the  irreversible  convection  losses 
through  the  side  walls  of  semiconductor  legs.  Higher  the  convection 
heat  transfer  coefficient,  more  heat  is  carried  away  without  being  con¬ 
verted  to  electricity.  A  certain  portion  of  the  heat,  entering  the  top 
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Fig.  3.  (a)  Contours  of  heat  transfer  and  (b)  i 
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tive  heat  transfer  boundary  condition  (h  =  15  W/m2K). 
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Fig.  4.  (a)  Contours  of  heat  transfer  and  (b)  contours  of  electric  potential  with  convective  heat  transfer  boundary  condition  (Ji  =  35  W/m2I(). 


surface  of  the  TEG,  is  leaving  through  the  side  surfaces  by  convection  as 
evidenced  by  the  heat  flux  lines  terminated  at  the  vertical  side  surfaces 
of  the  semiconductor  legs.  The  surface  of  a  semiconductor  leg  is  cooler 
than  the  core  due  to  the  heat  removal  by  convection  for  a  given  x  loca¬ 
tion.  The  heat  flux  lines  are  no  longer  parallel  in  both  p-type  and  n-type 
legs  and  a  non-uniform  distribution  of  heat  flux  line  is  observed.  Due  to 
the  coupled  thermo-electric  effect,  a  potential  difference  is  established 
in  the  TEG  which  drives  an  electric  current  through  the  system.  The 
electric  potential  remains  nearly  constant  at  a  given  location  of  distance 
x  when  h  =  0  for  a  particular  semiconductor  leg.  However,  non-linear 
electric  potential  distribution  is  observed  with  convective  heat  transfer 
boundary  condition. 

In  addition  to  the  numerical  simulation,  a  comparison  between 
analytical  result  and  numerical  simulation  is  presented  in  terms  of  the 
current  produced.  The  analytical  results  are  obtained  from  the  mathemat¬ 
ical  model  developed  by  Reddy  et  al.  [26],  The  comparison  is  presented  in 
Fig.  6(a)  where  the  produced  current  is  plotted  as  a  function  of 


temperature  difference  between  hot  and  cold  surfaces  of  the  TEG.  The 
electric  current  increases  as  temperature  difference  increases.  Fig.  6(a)  es¬ 
tablishes  a  good  agreement  between  the  current  result  and  result  avail¬ 
able  in  the  literature.  The  proposed  nanostructured  TEG's  thermal 
efficiency  is  demonstrated  in  Fig.  6(b)  where  thermal  efficiency  is  plotted 
as  a  function  of  temperature  difference  between  hot  and  cold  surfaces  of 
the  TEG  at  different  values  of  convection  heat  transfer  coefficients.  TEG 
has  highest  efficiency  with  larger  temperature  difference  and  adiabatic 
boundary  condition.  As  convection  heat  transfer  coefficient  increases, 
the  thermal  efficiency  drops  and  this  can  be  attributed  to  heat  loss 
shown  in  temperature  field  plots  from  Figs.  3(a),  4(a),  and  5(a). 

4.  Conclusion 

In  this  research  paper,  a  numerical  simulation  of  nanostructured  TEG 
is  carried  out.  The  nanostructured  TE  materials  have  low  thermal 
conductivity  and  higher  power  factor  (o ?cr)  which  improves  the 
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Fig.  5.  (a)  Contours  of  heat  transfer  and  (b)  i 
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!  boundary  condition  (h  =  50  W/m2K). 
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Fig.  6.  (a)  Comparison  of  current  production  using  numerical  and  analytical  techniques,  (b)  Thermal  efficiency  of  TEG  as  a  function  of  convection  heat  transfer  coefficient  and  temperature 
difference. 


performance  of  the  TEG.  Current  numerical  simulation  considers  the 
Seebeck  effect,  Peltier  effect,  Thomson  effect,  Fourier  heat  conduction, 
and  convection  heat  transfer  losses.  The  influences  of  hot  surface  tem¬ 
perature  and  convection  heat  transfer  coefficient  on  the  performance 
parameters  of  TEG  such  as  thermal  efficiency  and  electric  current  have 
been  studied.  Electric  current  generation  using  numerical  simulation 
and  analytical  simulation  shows  a  good  match.  An  increase  in  hot  sur¬ 
face  temperature  leads  to  increase  in  electric  current  generation  and 
eventually  the  thermal  efficiency.  Numerical  results  prove  that  the  pres¬ 
ence  of  irreversible  convection  heat  transfer  causes  a  large  amount  of 
heat  loss  thus  reducing  the  thermal  efficiency.  In  future,  more  detailed 
three-dimensional  numerical  simulation  of  TEG  will  be  carried  out  to 
observe  the  above  mentioned  effects  in  more  detail. 
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